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REINHARD ALKOFER 

FUNCTIONAL APPROACHES TO CONFINEMENT 
On the question of the time-Uke quark-gluon coupUng 

Jeff Greensite, when inviting to this round-table-dis- 
cussion, raised the question: What don't we know about 
Confinement? Although this question might seem sur- 
prising it is very well-posed. 

Let me follow for one moment the usually employed 
Gedankenexperiment of kicking a quark out of a hadron. 
As most of you know such reasoning leads to the de- 
velopment of a picture where a flux tube or string is 
formed which pulls either the quark back to the hadron 
or leads via string breaking and de-excitation to the 
production of (many) hadrons. 

As clever as these arguments are presented, typically 
an aspect is overlooked which, however, is only relevant 
in view of the available non-perturbative techniques. Lat- 
tice calculations and functional approaches rest almost 
exclusively on a formulation in Euclidean space-time. 
Processes which take place at time-like momenta are 
only describable with some additional effort. Given the 
fact that practically all physical processes involve time- 
like momenta this is of course a severe short-coming. 
E.g. taking your prefered quark propagator with "built-in 
confinement" might work nicely when considering static 
(or space-like) quantities but may fail miserably when 
calculating production processes. (NB: A corresponding 
example can be found in ref . [ 1 ] .) 



An understanding of confinement can probably only 
be considered complete if it includes knowledge about 
the coupling of gluons and quarks at time-like momenta. 
One can only appreciate the infeasibility of the required 
task when comparing to recent studies of the fully renor- 
malized quark-gluon vertex function at space-like mo- 
menta, see e.g. ref. [2] and references therein. First of 
all, one needs to notice that dynamical breaking of chiral 
symmetry also changes the character of the quark-gluon 
coupling: Not only vector-type chiral-symmetry preserv- 
ing couplings are allowed but also such type of couplings 
which are not invariant under chiral transformations, es- 
pecially scalar ones. Therefore, in order to discuss the 
issue of Lorentz scalar confinement we have to consider 
confinement and dynamical breaking of chiral symme- 
try in one unified treatment. Mentioning chiral symmetry 
one gets reminded that of course also the Ua(1) anomaly 
might very well be related to confinement, either directly 
(see e.g. [3]) or indirectly due to the fact that some topo- 
logically non-trivial field configurations are the "confin- 
ers" (see e.g. [4, 5, 6] and others). As the anomaly reflects 
itself in the tj' mass and in certain decays, again we will 
only truly understand whether it is correlated with con- 
finement when we gain some insight into the coupling of 
quarks and glue at time-like momenta. 

Therefore, either if we stay on purely theoretical play- 
grounds like investigating strong-coupling phenomena 
or we are really getting serious in calculating hadronic 
decays and/or production processes in terms of quarks 
and glue: An insight into their time-like properties and 
their respective mutual coupling is required. This task 



certainly demands to generalize and enlarge the available 
"tool box" of non-perturbative methods. 

A remark on suitable gauges 

There are not so many non-perturbative tools available 
to treat Quantum Field Theories at strong coupUng, and 
hereby especially the strong interaction as described by 
QCD. And there is basically only one method which does 
not require gauge fixing: Lattice Monte-Carlo calcula- 
tions. Corresponding calculations of e.g. hadron spectra 
are certainly very useful but without refering to the quark 
and glue substructure one is very restricted in learning 
about this substructure. As quark and gluon correlation 
functions are only non-vanishing if a gauge is fixed it 
can be very reasonable to fix a gauge. 

Acknowledging the reasons for fixing a gauge (which 
is a pure necessity in perturbation theory and functional 
approaches!) the following question arises: Are there 
gauges which are more suitable than others? This ques- 
tion has already been answered affirmatively by Hugo 
Reinhardt in this discussion, and I want to re-emphasize 
it again. 

The Landau gauge has proven to be especially suited 
for the investigation of Green functions, a guide to the 
literature can be e.g. inferred from the reference list of 
[7]. Even given the recent discussion on different types 
of solutions (which has been also very present at this 
Conference) one should not lose sight of how much 
knowledge has been gained in the last decade. 

Coulomb gauge is used in many investigations (see the 
references in [8, 9]) and allowed for some fundamental 
insight [10]. Therefore it is certainly worth investigating 
the properties of quarks and gluons in this gauge, espe- 
cially, as there is some recent progress after Coulomb 
gauge QCD has proven in the past to be much more com- 
plicated than expected, and this already at the perturba- 
tive level [11]. 

Last but not least, 1 would like to mention Maximally 
AbeUan gauge, and especially 1 like to point out that 
an infrared analysis of functional equations has been 
performed recently [12]. And, in addition to lattice 
calculations in this gauge (see e.g. [13] and references 
therein), also the functional equations are a subject of 
current investigations [14, 15]. This will allow further 
insight into the conjecture of Abelian dominance and, in 
the best of all cases, into the importance of chromomag- 
netic degrees of freedom as e.g. monopoles and center 
vortices. 

A comment on the infrared beliaviour of Landau 
gauge Green functions and lattice calculations 

As I already mentioned the debate on the infrared 
behaviour of QCD Green functions in the Landau gauge I 



would Uke to take the freedom to ask a question which is 
probably only understandable to the experts. (Therefore 
all non-experts might wish to skip this subsection.) 

It is usually claimed that lattice calculations (in four 
space-time dimensions) see only the so-called decou- 
phng solution (sometimes also called the "massive" so- 
lution [16]) of functional equations and not the scaUng 
solution. This is certainly true if one takes into account 
only the available data for the gluon propagator in the 
usual minimal Landau gauge. On the other hand, going 
to the extreme strong-coupling limit not only the infared 
exponents of the decoupling but also the exponents of 
the scaling solution can be extracted [17, 18, 19]. To my 
opinion such an extraction should be absolutely impossi- 
ble if the scaling solution is definitely absent and there- 
fore the lattice calculations should not contain any trace 
of it. But as both type of exponents are seen there is still 
something to be understood here, especially in the view 
of the fact that there are dependencies on the way the 
gluon field is discretized [17]. 

To summarize my contribution: Although the last 
years have seen tremendous progress in our understand- 
ing of the infrared behaviour of QCD almost all of the de- 
cisive questions cannot be definitely answered yet. How- 
ever, if one acknowledges also our progress in devel- 
oping new methods which are applicable in the strong- 
coupling domain I am convinced that in the not-so- 
distant future we will gain a basic understanding of the 
hardest problem of hadron physics, we will eventually 
understand confinement. 

DIVIITRI DIAKONOV 

HOW TO CHECK THATDYONS ARE AT WORK? 

There is multiple evidence from theory, lattice and phe- 
nomenology that dyons play an important role in en- 
acting confinement, and in inducing the deconfinement 
phase transition: 

• In ^ = 1 supersymmetric theory it is precisely 
dyons that shape the vacuum, in particular the 
gluino condensate, and it is an exact result there [20, 
21, 22] 

• Lattice studies show that zero fermion modes 
'jump' from one position in space to another as 
one varies fermion boundary conditions - precisely 
as one would expect from dyons carrying zero 
modes [23] 

• a semiclassical picture of the vacuum populated by 
dyons gives an appeahng explanation of all main 
features associated with confinement: the area be- 
havior of large Wilson loops, the asymptotic Un- 



ear rising potential only for nonzero A^-ality probes, 
the cancelation of gluons in the free energy, and 
a order deconfinement phase transition for all 
gauge groups except SU (2), regardless of whether 
the group has a nontrivial center or not [24, 25]. 

In a few words, if the Polyakov line is, on the average, 
not an element of the group center, as in the confinement 
phase, dyons appear as saddle points of the Yang-Mills 
(YM) partition function. Dyons are gluon field configura- 
tions with asymptotic Coulomb-like chromo-electric and 
-magnetic fields. The ensemble of many dyons is simi- 
lar to a multi-component plasma. In particular, the dual 
(magnetic) gluons obtain a Debye mass. This is the phys- 
ical reason for the exponential decay of the Polyakov 
lines correlations, i.e. of the hnear rising potential for 
heavy probe 'quarks' . The appearance of a mass for dual 
gluons is also responsible for the area behavior of large 
Wilson loops. A surface spanning the loop is a source for 
a soUton of finite thickness, 'made' of the dual fields. The 
action per area of this soliton is the string tension which 
coincides, at low temperatures T, with the string tension 
computed from the correlation function of the Polyakov 
lines [24]. 

The ensemble of dyons induce a nonperturbative en- 
ergy that is a function of the Polyakov loop eigenvalues. 
The function is such that its minimum, for any gauge 
group, is at a universal set of eigenvalues related to 
the Weyl vector, the half-sum of positive roots of the 
Lie algebra for the gauge group [25]. For most gauge 
groups it impUes that the trace of the Polyakov loop is 
zero in the lowest dimensional representations (however 
there are subtleties related to the fluctuations about the 
minimum). This dy on-induced nonperturbative potential 
energy competes with the well-known perturbative 
potential energy, also a function of the Polyakov loop 
eigenvalues. The perturbative energy scales as 7"* with 
respect to the nonperturbative one. Therefore, at some 
critical Tc it prevails, and that is the mechanism for the 
deconfinement phase transition. It happens irrespectively 
of what is the center of the gauge group. 

There are presently several qualitative pictures of con- 
finement being discussed. Apart from dyons, these are 
AbeUan monopoles and center vortices, see e.g. [26]. It 
may be that all three pictures can be, in a sense, recon- 
ciled. Dyons may be the real physical objects that reveal 
themselves as AbeUan monopoles seen on the lattice in 
the Abelian gauges, such as the maximal Abelian gauge, 
and also as center vortices observed in the center gauges, 
such as the maximal center gauge. 

Indeed, in order to assemble many dyons together they 
need to have the same asymptotic field A4 at spatial in- 
finity. That necessarily requires that dyons are in the 
'stringy' gauge where a singular Dirac string is stick- 



ing out from each dyon. The Dirac strings are gauge ar- 
tifacts: the action density is finite there. However, they 
do carry a quantized Abehan magnetic flux. When, in a 
lattice simulation, one uses any variant of the Abelian 
gauge, one identifies lattice magnetic monopoles as the 
sources of that flux. The exact position of Abelian lattice 
monopoles varies somewhat as one varies gauge fixing - 
in accordance with the fact that the direction of the Dirac 
string sticking from a dyon is subject to a gauge choice. 
However, lattice magnetic monopoles may well be a re- 
flection of real physical objects, the dyons. It would be 
interesting to check it directly. 

Furthermore, if one further restricts the gauge to a cen- 
ter gauge, center vortices are revealed [26]. They can be 
understood as the Dirac strings connecting dyons. There 
are gauges where a dyon has a Dirac string entering it, 
and another leaving it. See a recent study in Ref. [27] of 
the relation between dyons and vortices. 

1 should also mention a talk by Langfeld and llgen- 
fritz at this conference [28], who "cooled" lattice config- 
urations keeping Polyakov loops fixed. Usually smear- 
ing the configurations by cooling kills confinement but in 
this study it is preserved. The interesting observation is 
that the "cooled" configurations preserving confinement 
are mainly (anti)self-dual fields. Instantons and dyons are 
(anti)self-dual. 

Concerning instantons, the quantum Coulomb interac- 
tions of dyons are such that they tend to glue up into 
electric- and magnetic-neutral clusters which at low tem- 
peratures are hardly distinguishable from instantons [6]. 
The difference with the old 'instanton hquid' model [29] 
is that (0 the Polyakov line is now nontrivial, (»') the in- 
tegration measure over collective coordinates is invariant 
under permutation of dyons 'belonging' to different in- 
stantons, and allows instantons to overlap. These circum- 
stances are critical for obtaining confinement that was 
absent in former old instanton models. 

Finally, I would like to point out a (simple) lattice mea- 
surement that may help to understand the nature of the 
YM vacuum in general, and to demonstrate (or refute) 
the importance of dyons, in particular. I suggest to mea- 
sure the effective action for the gauge-invariant eigenval- 
ues of the Polyakov line. The definition is given in Eq.(3) 
of Ref. [25]. In lattice setting, one puts all time links 
to be unity matrices (the A4 = gauge), but make the 
spatial links periodic up to a gauge transformation with 
the matrix L(x) being the Polyakov line. Without loss of 
generality one can take it to be diagonal, as in Eq.(l) of 
Ref. [25]. One then simulates the ensemble of configura- 
tions with fixed 'eigenphases' 0(x) as boundary condi- 
tions. In particular, one can take (p to be x-independent. 
The partition function or the free energy itself is not cal- 
culable by Monte Carlo methods but one can find the av- 
erage plaquette and then integrate it over J3 or tempera- 
ture to obtain the free energy. This wiU be the effective 



potential as function of (|). At large T it is the well-known 
perturbative potential energy as function of ^. It is inter- 
esting to see how it looks like below and above Tc for 
different groups. 

If, as we assume, dyons are of relevance, the minimum 
of the effective potential will be at a specific point <^ 
proportional to the Weyl vector p, for any gauge group, 
seeEq.(5)ofRef. [25]. 

It should be noted that if one just studies the distri- 
bution of the Polyakov line 'eigenphases', it will be in 
any case dominated by the Haar measure weight that 
governs the local ultraviolet quantum fluctuations. In 
order to see the smooth potential energy as function of 
^, one really needs to consider the case of a constant or 
slowly varying Polyakov lines using, for example, the 
setting described above. 

JAN PAWLOWSKI 

CONFINEMENT AND THE PHASE 
DIAGRAM OF QCD 

An even partial answer to Jeff's question "What don't 
we know about confinement?" certainly gives rise to 
new ones, as is the case for all good questions. In this 
Uttle coiimient I would rather like to dwell on these 
new questions, as well as on weakened versions of the 
original question. 

In my opinion the question "What don't we know 
about confinement?" is closely related to the question 
"How do we detect confinement?". In Yang-Mills the- 
ory the latter question seems to be simple to answer: 
watch out for the Unear potential between static quark 
sources! Such a linear potential can be extracted from 
e.g. the area law for Wilson loops for large areas, or, at fi- 
nite temperature, from the large distance behavoir of the 
correlation function of a Polyakov loop and an adjoint 
Polyakov loop, just to name two of the many possibili- 
ties. This is fine as long as one does not add the ques- 
tion of the confinement mechanism. I think it is precisely 
here where the problems start. The above-mentioned op- 
erators are inherently non-local, and all attempts to de- 
scribe the confinement mechanism are based on at least 
semi-local degrees of freedom, mostly effective theories 
of topological defects. Loosely speaking one could say 
that the more local the description gets the less convinc- 
ing it is. In any case a quasi-local description implies a 
specific parametrization of the theory (a choice of a coor- 
dinate system), at least if quantum fluctuations are taken 
into account. Evidently such a choice is closely related to 
a gauge fixing procedure, and it is a hard-learned lesson 
in the past that the quest for the mechanism of confine- 



ment goes hand in hand with explicit or implicit gauge 
fixings as well as the question of their interrelation. 

This leaves us with the question of how to detect con- 
finement, in particular in formulations of QCD which are 
explicitly gauge-fixed. Again there have been a plethora 
of suggestions, for example based on the construction of 
the heavy quark potential, or on signatures of the absence 
of colored states in the physical Hilbert space. More re- 
cently much progress has also been made in directly con- 
structing observables that are sensitive to the confining 
properties of QCD. 

Still, in my opinion the above question is not fully 
clarified yet, and this fact is hampering the progress in 
our understanding of confinement in QCD with dynam- 
ical quarks. Due to string breaking the above definitions 
do not work anymore, and a clear signature of confine- 
ment is missing. This problem gets even more severe if 
we switch-on chemical potential/finite density. Then the 
question would be first "What does confinement mean 
in a dense medium of quarks/hadrons?" and then fol- 
lowed by "How do we detect it?". Again these apparently 
simple questions are very difficult to answer, in particu- 
lar on the basis of the other unresolved questions posed 
above, starting with Jeff's original one. Of course one is 
attempted to discard the latter questions as superfluous 
as we expect in this case the matter properties of QCD 
to dominate the physics. However, there seems to be a 
rather tight relation between chiral symmetry breaking 
and confinement (however we detect it) that even per- 
sists at finite density. In conclusion we will not under- 
stand this relation if we do not learn of how to answer 
the above questions about the glue sector in detail. 

Let me finish this little comment with toning down 
Jeff's original question to "What do we want to know 
about confinement?" and by collecting a couple of the 
answers given here: 

• What is the confinement mechanism? 

• How do we detect confinement? 

• How does confinement relate to chiral symmetry 
breaking? 

• What does confinement mean in a dense medium of 
quarks/hadrons? 

HUGO REINHARDT 

CONFINEMENT - WHAT WE DO KNOW, 
AND WHAT WE DO NOT KNOW 

Confinement is one of the millenium problems formu- 
lated by the Clay Mathematical Institute. To receive the 

corresponding prize requires to prove besides confine- 
ment, i.e. the absence of quarks and gluons from the 



physical spectrum, the existence of a mass gap and spon- 
taneous breaking of chiral symmetry. In the physics com- 
munity the generally accepted confinement criterium is 
that the temporal Wilson loop (order parameter) devel- 
ops an area law 

(iy(C)) = ^trexp(^-^A^^ -exp(-aE(C)) (1) 

or equivalently the spatial 't Hooft loop [30] (disorder 
parameter) with continuum representation [31] 

{V{C)) = (|exp (^-J ^(C)n) ^ ~ exp(-K:L(C)) 

(2) 

develops a perimeter law. Here, n(x) = 5/i5A{x) is the 
momentum operator of the gauge field and ^ (C) is the 
gauge field of a center vortex whose magnetic flux is 
located at the loop C. 

What we do know 

The 't Hooft loop operator in eq. (2) is a center vor- 
tex generator, which indicates that confinement must be 
related to the center of the gauge group. Indeed, from lat- 
tice calculations we know that in an SU (N) gauge theory 
an asymptotic string tension is obtained in all representa- 
tions with odd A^— ality, i.e. which transform non-trivially 
under the center of the gauge group [4]. Though we do 
not yet fully understand confinement, several pictures of 
confinement have emerged, which are supported by lat- 
tice results, see table 1 . 

The first two pictures rely on the condensation of mag- 
netic monopoles and center vortices, respectively, while 
the Gribov (-Zwanziger) picture requires an infrared di- 
verging ghost form factor d defined by 

{{-Dd)-')=d/{-A). (3) 

Here —Dd is the Faddeev-Popov operator in Landau or 
Coulomb gauge with D being the covariant derivative in 
the adjoint representation. 

Empirically it was found on the lattice that the dif- 
ferent pictures of confinement are Unked to each other: 
Magnetic monopoles are located on center vortices, [4], 
and both configurations are on the Gribov horizon in 
both Coulomb and Landau gauge [32]. Removing cen- 
ter vortices removes the magnetic monopoles and makes 
the ghost form factor infrared finite [33]. Furthermore, in 
the absence of magnetic monopoles, center vortices are 
topologically trivial in the sense that they have vanishing 
Pontryagin index [34]. Topologically trivial field config- 
urations cannot account for the topological susceptibiUty 
and spontaneous breaking of chiral symmetry. 

We have lattice evidence that chiral symmetry break- 
ing is lost when the "confiners" Uke center vortices and 



TABLE 1. Survey of confinement scenarios 



scenario 


dom. field config. 


gauge 


dual supercond. 


mag. monopoles 


Max. Abelian 


vortex cond. 


center vortices 


Max. Center 


Gribov- 


Gribov horizon 


Landau/Coulomb 


Zwanziger 







magnetic monopoles are removed from the Yang-Mills 
ensemble [4]. By the Banks-Casher relation [35] a non- 
vanishing quark condensate (the order parameter of chi- 
ral symmetry breaking) requires a non-zero quark level 
density at zero virtuality. As shown in ref. [36] the quark 
spectrum develops a gap at zero virtuality when center 
vortices are removed from the Yang-Mills ensemble. 

We also know that confinement is realized only in the 
low temperature and low matter density phase and a de- 
confinement phase transition occurs when the tempera- 
ture or matter density exceeds certain critical values. Fur- 
thermore, in the center vortex picture (and similarly in 
the monopole picture) the deconfinement phase transi- 
tion is a depercolation transition [37]. 

What we do not know 

So far, we do not have a rigorous proof of confinement 
in Yang-Mills theory in the sense of the Clay Mathemat- 
ical Institute. Furthermore, we have no gauge invariant 
confinement picture. What is common to all the confine- 
ment pictures hsted in table 1 is that they rely on spe- 
cific gauges. Of course, confinement is a physical phe- 
nomenon, which does not depend on a specific gauge 
choice. However, the confinement mechanism may be 
easier to reveal in one gauge than in another. There is 
nothing wrong in choosing a convenient gauge (choos- 
ing a gauge is nothing more than choosing specific "co- 
ordinates") but one has to show that the confinement 
picture found in one gauge persist when the gauge is 
changed. The ultimate confinement scenario should be 
formulated in a gauge invariant way. It is clear that the 
Gribov-Zwanziger scenario intrinsically requires a gauge 
fixing. Furthermore, the Gribov-Zwanziger mechanism 
seems to be realized in Coulomb gauge but not in Landau 
gauge. We also do not know the nature of the deconfine- 
ment phase transition in Gribov's confinement scenario. 
Clearly, the Gribov horizon of Coulomb gauge is inde- 
pendent of the temperature. 

Finally, we do not know yet how spontaneous break- 
ing of chiral symmetry is accomplished in the various 
pictures of confinement. If confinement and chiral sym- 
metry breaking are triggered by the same mechanism 
then the restoration of chiral symmetry should precisely 
occur at the deconfinement phase transition. This is 
confirmed on the lattice [38], [39]. 



Results of the Hamilton approach to Yang-Mills 
theory in Coulomb gauge 

During the last years our group has worked on Yang- 
Mills in Coulomb gauge, both in the continuum and 
on the lattice and the obtained results strongly sup- 
port Gribov's scenario of confinement. In the continuum 
an approximate variational solution of the Yang-Mills 
Schrodinger equation was carried out with the following 
results: 

1 . The gluon energy co (^) is infrared divergent and can 
be nicely fitted by Gribov's formula 

(0{k) = ^JP+W^/P , (4) 

where M w S60MeV [40], [41]. The same result 
is also obtained on the lattice [42]. Hence, gluons 
are absent from the physical spectrum, which is a 
manifestation of gluon confinement. 

2. The 't Hooft loop is found to obey a perimeter law 
[43]. 

3. A Unear rising non-AbeUan Coulomb potential was 
found [44], which is necessary but not sufficient for 
the area law in the Wilson loop. 

4. The inverse ghost form factor {k) was shown to 
represent the dielectric function of the Yang-MiUs 
vacuum [45]. As a consequence, the horizon con- 
dition d~^{k = 0) = implies that the Yang-Mills 
vacuum is a dual superconductor, which relates Gri- 
bov's confinement scenario with the monopole con- 
densation picture. 

The results of the variational approach were not only 
confirmed by lattice calculations but also by alternative 
approaches to continuum Yang-Mills theory like Dyson- 
Schwinger equations [46] and renormaUzation group 
flow equations [47]. 

VALENTIN ZAKHAROV 

SCALAR FIELDS AND CONFINEMENT 

Scalar, or Higgs fields are an indispensable part of the 
Standard Model. Another scalar field, the inflaton, is a 
common device used to explain cosmological inflation. 
The theory of conflnement is in the same sequence: it 
starts with an analogy to superconductivity and intro- 
duces a scalar field (pQ^ ^ la (relativistic version of) 
Ginzburg-Landau theory. In particular, for the tension of 
the string connecting heavy quarks one would expect: 

Cfstring ~ < <1>GL >^ , (5) 



In each case we deal with field theories which are, at 
first sight, well defined. And nevertheless it is common 
knowledge that in fact such theories are inconsistent 
because of the quadratic divergences. Thus, we have to 
expect the field theories to be replaced by something else. 

It is only in case of confinement, however, that there 
exists an independent source of information on the nature 
of the scalar field, that is lattice simulations. And, indeed, 
the interpretation of the data (Chernodub and Zakharov, 
2003) looks unusual: 

< |<P™jn| > ~ ^QCD^ (6) 
(< <Pmagn >)^ ~ Ag^^a , 

where (pmagn is the magnetically charged scalar responsi- 
ble for confinement, as seen on the lattice, a is the lattice 
spacing, a in the continuum Umit. Thus, in contrast 
to (5), the vacuum expectation value vanishes in fact and, 
nevertheless, suffices to ensure confinement at any a^Q\ 

From dimensional considerations alone it is obvious 
that the Eqs (6) are indicating to strings, or 2d field the- 
ory. Moreover, it is also seen in the simulations that the 
scalar-particle trajectories cover densely 2d surfaces, or 
strings (defined independently). Thus, the UV quadratic 
divergence seems to be avoided by invoking strings. (It is 
remarkable that it is known since long (Atick and Witten, 
1988) that strings are much 'softer' in terms of effective 
degrees of freedom than field theories. In that case, two 
powers of temperature were missing in the limit of high 
temperatures. At T = 0, two powers of the UV cut off are 
absent from < >.) 

In the continuum language, the most recent candidate 
for a 'soft' scalar field is the Horava-Lifshitz scalar, 
(Horava, 2009). It is assumed that at short distances 
theory becomes time-space asymmetric, L ~ (<3;0)^ — 
(A^)^. Furthermore, one argues that the effective field 
theory is in fact two-dimensional at short distances. In 
this sense, the solution for the scalar is similar to what 
we have just discussed with reference to the lattice. 

The lattice data refer to the Euclidean case and there 
can be no time-space asymmetry. While the Horava- 
Lifshitz scalar refers to the Minkowski space. The match- 
ing can be discussed only in terms of matrix ele- 
ments, like < >, which can be safely continued to 
Minkowski space. If the matching observed is meaning- 
ful indeed, then there arises an intriguing possibility that 
confinement effects in Minkowski space and at short dis- 
tances do not observe Lorentz invariance, whatever it 
might mean. 



where ^ is a scalar field. 



DAN ZWANZIGER 

CAN THE FREE ENERGY WITH SOURCES 
BE MEASURED ON THE LATTICE? 

Jeff Greensite has posed to the panelists the question 
"From your own perspective, what would be the most 
useful piece of information that you could obtain (from 
either lattice experiments, real experiments, analytical 
methods, or maybe from God) that would advance your 
own program, and/or our general understanding of con- 
finement?" I will suggest a lattice measurement of a per- 
haps novel type, and explain why it appears useful. Be- 
fore doing so, I frankly confess that I don' t know whether 
the suggested measurement is feasible in practice. 

Recent numerical studies on large lattices of the gluon 
propagator in momentum space D{k), reviewed recently 
in [48], yield finite values for D{0), in apparent dis- 
agreement with the theoretical expectation that D{0) = 0, 
originally obtained by Gribov [40], and argued in [49]. 
Upon reviewing the argument [49] which leads to this 
result, one hypothesis stands out which should perhaps 
be dropped in view of the apparent disagreement. This is 
the hypothesis that the free energy W{J) in the presence 
of sources J is analytic in J. This is an important point 
because a non-analyticity in the free energy is character- 
istic of a change of phase. 

The free energy W{J) enters the picture because it is 
the generating function of the cormected gluon correla- 
tors. In particular the gluon propagator Dj^^, is a second 
derivative ofW{J) at 7 = 0, 



dJrdJv 



7=0 



(7) 



Here a condensed index notation is used, where Jx repre- 
sents (x), ^ is a Lorentz index, and bis a color index, 
and we write 



[J, A) = I^A, = / J^x 7^(x)A^(x). 



(8) 



The free energy W{J) in the presence of sources J is 
given by 

expW{J) = (exp(y,A)) 

= / dAp{A)exp{J,A). (9) 

For simplicity we have written this in continuum theory, 
but we have in mind a numerical study on the lattice, and 
the lattice analog of each formula should be obvious. The 
integral over A is effected in Landau gauge d^A^ = 0. 
The domain of integration is restricted to the Gribov 
region Q., a region in A-space where the Faddeev-Popov 
operator is non-negative, M(A) = —dnD^{A) > 0. In 



practice, the positive, normalized, Euclidean probability 
density p(A) will depend on the numerical algorithm 
which is used for gauge fixing. 

[To give a concrete example, one may imagine an 
idealized gauge fixing - not attainable in practice - to the 
absolute minimum of the minimizing functional /^(g) = 
I |^A| p on each gauge orbit. In this case, in the absence of 
quarks, W{J) is given by 

expW(/)= / t/Ap(A)exp(7,A), (10) 
Ja 

where the fundamental modular region A C is a subset 
of the Gribov region, p(A) is given by the Faddeev- 
Popov formula, 

p{A)=N5id •A)det(M(A))exp[-5yM(A)], (11) 

and Sym{A) is the Yang-Mills action.] 

The gluon n-point function is the n-th derivative of 
W{J) at 7 = 0. To elucidate the analytic structure ofW{J) 
it would be useful to measure W (J) itself for finite values 
of J, instead of only of few of its derivatives at7 = 0. For 
example, one could take 



4{x)=H^^cos{k-x), 



(12) 



and measure 



sxpW{k,H) = { exp[£//^cos(A;-x)A^(x)] ), (13) 

X 

for a set of values of the parameters A:^ and //^. By 
analogy with spin models, H^ may be interpreted as the 
strength of an external magnetic field, with a color index 
b, which is modulated by a plane wave cos{k ■ x). The 
measurement is done using the standard gauge-invariant 
Wilson action followed by gauge fixing to the Landau 
gauge by a minimization algorithm. 

In [49] it was shown that the free energy per unit 
(Euclidean) volume satisfies the bound 



W{k,H) 
V 



<c\k\\H\ 



(14) 



where c is a certain constant. It would be of interest to 
investigate numerically whether this bound is saturated 
asymptotically at low k. If so, it would be a remark- 
able non-analyticity because the bound is linear in \H\, 
whereas normally one expects W{k,H) to be a power se- 
ries in H, with leading term of order H^. However such 
analytic behavior leads to the above-mentioned disagree- 
ment with the lattice data. 

It may be that importance sampling based on p(A) 
gives poor results for (exp(7,A)) because of the expo- 
nential character of the observable exp(7,A), although it 



is in fact real and positive. If so, one could alternatively 
try to measure A{x) in the presence of finite source J, 

using importance sampling with the positive weight 

p^(A)=p(A)exp(7,A). (16) 

This has the difficulty that the importance sampling must 
be done after gauge fixing, using the Faddeev-Popov 
weight 

p(A) = 5{d •A)det(M(A))exp[-5YM(A)]. (17) 
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